Specific heat of Bao esKo 32Fe2As2: evidence for multiband strong-coupling superconductivity 
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The specific heat of high-purity Bao.68Ko.32Fe2As2 single crystals with the highest reported superconducting 
Tc = 38.5 K was studied. The electronic specific heat, Cp, below Tc shows two gap features, with Ai ~ 11 
meV and A2 ~ 3.5 meV obtained from an Q-model analysis. The reduced gap value, 2A'^'^^ /kBT^ ~ 6.6, the 
magnitude of the specific heat jump, ACp(Tc)/Tc, and its slope below Tc exhibit strong-coupling character. 
We also show that an Eliashberg model with two hole and two electron bands gives the correct values of Tc, the 
superconducting gaps, and the temperature dependence of the free-energy difference. 



The newly discovered iron pnictide superconductors 
present an unusual case of multiband superconductivity in 
the vicinity of a magnetic instability. In view of the weak 
electron-phonon coupling in this class of superconductors \ T\, 
magnetic excitations are the most promising candidates for 
the pairing boson. The analysis of experimental data and the 
theoretical modeling of these systems are, however, consid- 
erably complicated by their complex electronic structure that 
involves at least four energy bands crossing the Fermi surface. 
Despite intensive research, there is hence no consensus on the 
gap symmetry and on the nature of the pairing interaction. 

The thermodynamic properties are a key source of surface- 
insensitive information on the electronic interactions and on 
the structure of the superconducting gap function. For a su- 
perconducting order parameter with alternating sign, which 
appears to be present in at least some of the pnictides, the 
pairing interactions are strongly influenced by impurity scat- 
tering, so that sample quality and phase purity issues are cru- 
cial for calorimetric experiments. Presumably because of such 
issues, the current specific heat evidence of doped BaFe2As2, 
a pnictide family for which sizable single crystals are avail- 
able, is still ambiguous. Recently two gaps were reported in 
electron-doped Ba(Feo.925Coo.o75)2As2 , supporting the re- 
sults of nuclear magnetic resonance (NMR) and muon spin 
relaxation (/iSR) studies |2|. In hole-doped Bai_xKxFe2As2, 
on the other hand, prior specific heat experiments have shown 
only one gap with values of 6 meV in a single crystal with Tc 
= 36.5 K [3J and varying from 5.9 meV to 6.7 meV for 0.3 
< a; < 0.6 (36.0 K < i; < 37.3 K) in polycrystafline speci- 
mens 14], in contrast to the observation of two nearly isotropic 
gaps in angle-resolved photoemission (ARPES) experiments 
(|5|. Alas, all of the specific heat data thus far reported suffer 
from a residual low-temperature non-superconducting elec- 
tronic contribution and show Schottky anomalies. In electron- 
doped BaFe2As2, fliis complication is likely due to structural 
defects caused by Co substitution directly in the supercon- 
ducting FeAs planes. In hole-doped BaFc2As2, the substi- 
tution of Ba by K does not disturb the structure of the FeAs 
planes, but gives rise to electronic phase separation in the un- 
derdoped regime ||7] [S) . 

In this Letter we report a comprehensive study of the 
specific heat, Cp, of an optimally hole-doped, high-purity 



Bao.68Ko.32Fe2As2 single crystal with Tc = 38.5 K (onset), 
the highest transition temperature reported to-date for this 
family of superconductors. The superconducting specific heat 
indicates two energy gaps with magnitudes 2Ai — 6.6 ksTc 
and 2A2 = 2.2 ksTc, in agreement with flie ARPES, ^SR, 
and NMR results f5"6l. We also calculated the specific heat 
from the Eliashberg spectral function in a four-band spin- 
fluctuation model. From the solution of the Eliashberg equa- 
tions we obtained the value of Tc, the superconducting gaps, 
and the temperature dependence of the free-energy difference 
in good agreement with the experimental data. The averaged 
electron-boson coupling constant. A"" w 1.9, resulting from 
this analysis and consistent with the magnitude of the specific 
heat jump and its first derivative below Tc, is distinctly larger 
than the estimated electron-phonon coupling 1 1 1 . The strong 
coupling suggests the interaction with low-energy (< 50 me V) 
interband spin-fluctuations as the pairing mechanism. 

The superconducting Bao.68Ko.32Fe2As2 and reference 
Ba(Feo.88Mno.i2)2As2 single crystals were grown from self- 
flux in zirconia crucibles sealed in quartz ampoules under ar- 
gon atmosphere, as described earlier 1(91. Their chemical com- 
positions were determined by energy-dispersive X-ray spec- 
trometry. The specific heat was measured in the temperature 
range between 2 and 200 K with a Physical Properties Mea- 
surement System (Quantum Design) using the thermal relax- 
ation technique. The samples were mounted on a standard 
puck using Apiezon N grease. The measurements showed 
no difference between several cleaved pieces from the same 
batches. Bulk superconductivity in Bao.68Ko.32Fe2As2 was 
confirmed by magnetic susceptibility, dc and infrared conduc- 
tivity measurements, whereas Ba(Feo.88Mno.i2)2As2 did not 
show any signatures of superconducting or magnetic phase 
transitions. In the following, we present specific-heat data 
for a superconducting Bao.68Ko.32Fe2As2 crystal of weight 
13.6 mg and for a Ba(Feo.88Mno.i2)2As2 reference crystal of 
weight 9.9 mg. 

Figure la displays unprocessed specific heat data on the 
Bao.68Ko.32Fe2As2 and Ba(Feo.88Mno.i2)2As2 single crys- 
tals. Cp{T)/T of Bao.68Ko.32Fe2As2 shows no signs of low- 
temperature upturns or Schottky anomalies at any magnetic 
field. In the low-temperature limit (inset in Fig. la), the data 
can be fitted by C(T)/T = 7(0) + I3T^, where 7(0)r rep- 
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FIG. 1: (a) Temperature dependence of the specific heat Cp/T 
of Bao.68Ko.32Fe2As2 (squares) and Ba(Feo.88Mno.i2)2As2 (cir- 
cles). The inset shows a plot of Cp/T vs. at low T for both 
samples. The lines represent the best fit to Cp (T) /T = 7(0) + /JT^. 
(b,c) Temperature dependence of Cp of Bao.68Ko.32Fe2 As2 near 
measured at different external magnetic fields applied perpendicular 
to the ab-plane. 



resents the residual zero-temperature electronic specific heat 
and the second term the lattice specific heat. The fitting pa- 
rameters for zero field are 7(0) = 1-2(2) mJ/mol and 
/3=0.496(1) mJ/mol K". For magnetic field H = 9 T we used 
the same /3 and obtained 7(0) ~ 6.4(2) mJ/mol K^. These 
residual specific heat values are the lowest reported for FeAs- 
based superconductors so far, indicating the superior quality 
and high purity of our samples. The ratio of the residual elec- 
tronic specific heat to its normal-state counterpart (7^, see 
below) yields an estimate of the non-superconducting phase 
fraction 7(0)/7f « 1.2/50 = 2.4 %. 

Figure lb highlights the superconductivity-induced jump 
anomaly of Cp, which exhibits the largest magnitude (125 
mJ/mol K^) and smallest width (< 0.4 K) reported to-date |2 - 
m [To] [n] . An external magnetic field applied perpendicular 
to the a6-plane gradually suppresses the jump and reduces the 
transition temperature by 1 .4 K for H — QT. The transition 
temperature Tc determined as the maximum of the derivative 
of Cp (T) decreases linearly with magnetic field with a slope 
of 5Hc2{T)I5T\t= - 6.7 T/K (Fig. 2c), that is consistent 
with the value reported in Ref. (TO) This is a factor of two 
larger than the value obtained from resistivity measurements 
on samples from the same batch lfT2l . This difference may be 
due to flux-flow effects. 
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FIG. 2: Superconductivity-induced specific-heat difference of 
Bao.68Ko.32Fe2As2. The black solid curve is the result of a fit 
according to the two-gap a-model. The dashed and dash-dotted 
curves represent the partial contributions of the two bands. The light 
gray line is the electronic specific heat calculated from the four-band 
Eliashberg model. Inset: Temperature dependence of the thermody- 
namic critical field (see text). 

In order to reliably extract the electronic contribution, 
Cei{T), from the total measured specific heat, the contri- 
bution of lattice excitations, Ciatt{T), has to be accurately 
determined. The lattice specific heat of the parent com- 
pound BaFe2As2 can not be accurately obtained because of 
the magnetic and structural phase transitions at ^ 140 K. 
We found that substitution of 12% Fe by Mn suppresses the 
spin-density-wave state and does not induce superconductiv- 
ity. The specific heat of Ba(Feo.88Mno.i2)2As2 (red cir- 
cles in Fig. la) was found to be independent of magnetic 
field. The electronic and lattice specific heat terms extracted 
from a fit to low-temperature data on this compound (inset in 
Fig. la) are 7(0) = 7^,, = 14.9(2) mJ/mol and Pmu = 
0.420(4) mJ/mol K^, respectively. The lattice specific heat of 
Ba(Feo.88Mno.i2)2As2, C^^^l(T), was obtained by subtract- 
ing 7M„r from its total specific heat, C^^l^T). The esti- 
mated Debye temperature is 306 K, somewhat higher than the 
corresponding value for superconducting Bao.68Ko.32Fe2As2 
(277 K). Since the difference of the lattice parameters of both 
compounds does not exceed 1.5%, we assume that the phonon 
contributions to their specific heat and to the entropy obey a 
law of corresponding states. The normal-state specific heat of 
Bao.68Ko.32Fe2As2 can then be obtained from the commonly 
used corresponding states approximation [|i3|: 

CL{T) = Q^«(T)+7f -T = A-Ctlrt{B-T)+^f,-T (1) 

where 7^ is the Sommerfeld constant in the normal state of 
Bao.68Ko.32Fe2As2, and A and B are close to unity. From 
a least-squares fit of our data to Eq. (1) within the temper- 
ature range 40 - 150 K (more than 300 data points) under 
the constraint of entropy conservation, we obtained 7^= 50 
mJ/mol K^, with A = 0.95 and B = 1.03. Figure 2 shows the 
difference, ACei{T)/T, between the measured specific heat 
Cp{T)/T of Bao.68Ko.32Fe2As2 (Fig. la) and its normal- 
state counterpart C^i{T)/T estimated from Eq. (1). The 
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same ACe/ (T) /T (within the symbol size in Fig. 2) has also 
been obtained by representing the lattice contribution on a ba- 
sis of Einstein modes [14|. The value of the Sommerfeld con- 
stant 7^^= 50 mJ/mol is in between those reported for the 
pristine end compounds BaFe2As2 (6.1 mJ/molK^ lilSll ) and 
KFe2As2 (69.1 mJ/molK^ USD- This value exceeds the one 
calculated from the band structure (joft ~ 10.1 mJ/molK^ 
for Bao.6Ko.4Fe2As2 [17 1) by almost a factor of five, while 
it is consistently inferred from all specific heat experiments 
reported to date ISlffllTTI. 

The thermodynamic critical field, Hc{T), can be obtained 
from the free-energy difference between the normal and su- 
perconducting states: 

^F{T) - ^iQVmHl{T)/2 - A[/(T) - T^S{T) (2) 

where Vm = 6.05 x 10^^ m'^/mol is the molar volume 
determined from our X-ray diffraction measurements on 
Bao.68Ko.32Fe2As2. The temperature dependence of the ther- 
modynamic critical field obtained by numerical integration 
of our data is displayed in the inset of Fig. 2. By fit- 
ting the thermodynamic critical field expressed as Hc{T) = 

Hco (l - (T/Tcf'j to our data between 2 and 38.5 K we ob- 
tained i.ioHco = 0.85 T, a value comparable to those of the 
cuprate superconductors IfTSl . 

In contrast to 7^, the superconductivity-induced elec- 
tronic specific heat is very sensitive to the sample qual- 
ity and phase purity. We now show that if the impu- 
rity scattering is minimized, ACei(T) reveals the intrin- 
sic multigap strong-coupling nature of the superconduct- 
ing state. Two quantities are commonly used to assess 
the pairing strength, namely the reduced jump anomaly, 
ACe;(rc)/(7^Tc), and the normaUzed slope of ACe;(T) 
right below r„ 5(T,) = - d[ACeiiT)/ij^T,)]/dT |t,. 
Their magnitudes in Bao.68Ko.32Fe2As2, 2.5 and 11.1, re- 
spectively, are significantly larger than the corresponding BCS 
values of 1.43 and 3.77. In the framework of the Eliash- 
berg formalism, the coupling strength can be expressed in 
terms of the ratio hujin/kBTc, where aj;„ is the logarithmi- 
cally averaged frequency of bosons mediating the Cooper 
pairing |19|. In Figure 3, we place the reduced superconduct- 
ing jump and g-value of Bao.68Ko.32Fe2As2 on a universal 
plot 1 19 1 as a function of this dimensionless coupling param- 
eter Both values correspond to huJin/kBTc = 7.5, close to 
the prototypical strong-coupling superconductor Pbo.8Bio.2- 
This comparison indicates strong-coupling superconductivity 
inBao.68Ko.32Fe2As2. 

The same issue can be addressed by considering the mag- 
nitude of the superconducting energy gap. Clearly, the promi- 
nent knee in the ACei{T)/T data around T = 15 K (Fig. 2) 
cannot be described in terms of models with a single gap. We 
have hence fitted these data to the phenomenological multi- 
band a-model |20|, which assumes a BCS temperature de- 
pendence of the gaps and has been widely used to analyze heat 
capacity data on MgB2 121 1. The superconducting gap mag- 
nitudes at T = are introduced as adjustable parameters aj. 
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FIG. 3: Comparison of the dimensionless superconducting jump, en- 
ergy gap at T = 0, and slope of the specific heat near Tc calculated 
by the semiphenomenological expressions in Ref. 19 (lines) with 
our results for Bao.68Ko.32Fe2As2 (open circles), andforPbo.8Bio.2 
(triangles) and In (squares), after Ref. J9 The open and solid circles 
represent the maximum gap values obtained by the two-band a and 
four-band Eliashberg model, respectively. Shaded areas represent the 
uncertainty in uJin/Tc according to the error bars in our data points. 



defined according to Aj{T) = (a^ / aecs) Abcs> where 
Q!BCS = ABcs(0)/fc_B7c = 1.764 is the weak-coupling value 
of the gap ratio. Another set of adjustable parameters is the 
fractions of the total electronic density of states, 7j/7x, that 
each band contributes to the superconducting condensate. Fig- 
ure 2 shows the result of our multigap fit (black solid curve) 
with only two gaps survived ai = 3.3, a2 = 1.1, and 71 ^ 
72 0.5 • 7j^, which reproduces well ACei {T)/T below T^. 
One of the gaps, A2(0) = 3.5 meV, is somewhat smaller than 
Abcs(O), whereas the other one, Ai(0) = 11 meV, is much 
larger Our analysis thus yields thermodynamic evidence 
of two different superconducting gaps (or groups of gaps) 
with quite different absolute values in Bao.68Ko.32Fe2As2, in 
agreement with ARPES and /iSR results ||5l. Further, mark- 
ing the corresponding point, 2A'^^^ /kBTc = 6.6, on the uni- 
versal plot of Fig. 3 strongly supports our conclusion about 
the strong-coupling nature of superconductivity in this com- 
pound. 

Finally, we illustrate that the strong coupling scenario can 
be described by fully microscopic calculations in the frame- 
work of a spin-fluctuation model. Such models generically 
yield s-wave superconducting gaps with different signs on dif- 
ferent bands ("s± state") ll22l . It was shown in Refs. I23ll24l 
that strong-coupling effects lead to merging of the gaps in a 
simple two-band s± model. Therefore, a more realistic four- 
band model based on a band-structure with two hole bands 
and two electron bands crossing the Fermi surface needs to 
be considered |25 1. The main input into the Eliashberg equa- 
tions is the spectral function of the intermediate bosons. Fol- 
lowing Ref. |26] we took a spin-fluctuation coupling func- 
tion Bij{il) = XijBiVl) [inset of Fig. 4] with a linear VI- 
dependence at low frequencies, a maximum at hW^p^ = 18 
meV, and a fast decay at > hfl'gp'^, in qualitative agree- 
ment with recent experimental data on the normal-state dy- 
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FIG. 4: Free-energy difference between normal and superconducting 
states, experimentally obtained (open circles) and calculated from 
the four-band Eliashberg model (light gray line). The dashed, dash- 
dotted, and dotted lines represent the partial contributions of the in- 
dividual bands. Inset: Spin-fluctuation coupling function B{Q). 

namical spin susceptibility of Ba(Fei_xCox)2As2 ll27l l28l . 
Here Ay is the coupling constant for pairing of electrons in 
bands i and j. For our calculations we use the intraband cou- 
pling matrix elements Xu = 0.2 in order to take account of 
the weak electron-phonon contribution (T], for interband re- 
pulsion Ai2 — A34 = due to the symmetry of the wave 
functions, and A13 — A14 = —1.0, A23 — X24 = —0.2 
(Xji = XijN^''^{0)/Nj''^(0)). The corresponding densities 
of states are taken as N^{0) = 29 Ry \ iV^'(O) = 43 Ry^^ 
for the hole bands, and as iV|(0) = Nl{0) = 8.5 Ry^^ for 
the two equivalent electron bands (14]. The chosen param- 
eters, Xij{i 7^ j) and nI''''^{0), allow the best simultaneous 
fit to the experimental values of Tc, the superconducting gaps, 
and the temperature dependence of the free-energy difference, 
yielding = 38.5 K, and the following gap values: Aj* = 
-8.5 meV, = -3.6 meV for hole bands, and A;^ = A| 
= 9.2 meV for electronic ones. The effective coupling con- 
stant averaged over all bands, A°*' = J^ij ^i^'"^ / ^tou 
has a value of 1 .9, remarkably close to the coupling constant 
reported for Pbo.8Bio.2 and in agreement with the con- 
clusions of our phenomenological analysis presented in Fig. 
3. The consistency with the experiment tolerates some vari- 
ation of the parameters and remains satisfactory for MVg^ 
within the 10 - 20 meV energy range 1 14 1. A shift of MVgf^ to 
lower energy leads to a larger values of A"" with N^''^{0) ap- 
proaching the results of density functional theory (DFT) {VT\. 
This accounts for the strong renormalization of the Sommer- 
feld constant ^^Iidft ^ 5. 

The superconductivity-induced free-energy difference, 
AF{T), was calculated by using the expressions in Ref. l2T| 
The result is presented in Fig. 4, along with the partial band 
contributions to AF{T) which demonstrate that the super- 
conductivity in the second hole band has an induced origin. 
Fig. 4 also shows that the result of the model calculation is 
in fairly good agreement with the free-energy difference ob- 
tained by integrating the experimentally measured ACe/(T), 



while some deviations between the calculated and observed 
specific heat can be seen in Fig. 2. Some such deviations 
are expected because feedback effects of superconductivity on 
the bosonic spectral function, which lead to the formation of 
a temperature-dependent "resonant mode" in the spin fluctua- 
tion spectrum below Tc ll27l . have not been considered in the 
calculations. 

In summary, the specific heat of hole-doped 
Bao.68Ko.32Fe2As2 was measured on crystals of supe- 
rior quality and phase purity. Our a-model fit to the 
electronic part of the specific heat shows that there are two 
groups of superconducting gaps with magnitudes of 11 and 
3.5 meV. The magnitude and shape of the superconducting 
jump anomaly indicate strong coupling with intermediate 
bosons, in good agreement with an Eliashberg analysis of a 
spin-fluctuation model. 
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SUPPLEMENTARY MATERIAL 



• Evaluation of the normal-state specific heat at low temperature 

In order to reliably extract the electronic contribution, Cei{T), from the total measured specific heat, the contribution of 
lattice excitations, Ciatt{T), has to be accurately determined. The conventional method to evaluate the lattice specific heat of a 
superconductor is to recover the normal state by applying a magnetic field. This is not feasible for Bao.68Ko.32Fe2As2, because 
the requisite field range is experimentally inaccessible. In our work, by using the lattice specific heat of Ba(Feo.88Mno.i2)2 As2, 
C/a^M (T), as a reference, we evaluated C^^^j (T) of Bao.68Ko.32Fc2 As2 as A ■ Cl^H {B ■ T), with scaHng factors A and B of 0.95 
and 1.03, respectively. The green curve in Fig. SI represents the resulting normal-state specific heat of Bao.68Ko.32Fe2As2 from 
a least-squares fit of our data (open circles) to Eq.(Sl) within the temperature range 40 - 160 K under the constraint of entropy 
conservation. 

Alternatively, we also simulated Cl^^^{T) by representing phonon density of states (PDOS) F{lu) by six (5-functions Fi6{Ld — 
LJi) on a basis of Einstein frequencies in a geometrical ratio w^+i 
given by 



1.75 • Wi. The corresponding normal-state specific heat is 



1^ 



T 



1) 



1^ 



(SI) 



where R is the universal gas constant, Xi — hLUi/ksT, and the sum of the weights Fi is constrained to the number of atoms per 
unit cell. A justification of this trial PDOS is addressed, e.g., in Refs. S1-S3. The number of modes was chosen to be small 
enough to ensure the stability of the solution. The geometric ratio 1.75:1 allowed us to avoid spurious oscillations of the phonon 
DOS. The weights Fj , the starting Einstein frequency oji, and the Sommerfeld constant 7^ (« 49 mJ/molK^) were determined 
by a least-squares fit of the high temperature (40 K < T < 200 K) specific heat of Bao.68Ko.32Fe2 As2 (open circles in Fig. S2) 
to Eq.(Sl) under the constraint of entropy conservation. The results of the fit is given in Table SI. The (5-functions of the phonon 
DOS F{uj) are represented by a histogram in the inset of Fig. S2. The frequencies fuvi are ranging from 2 to 32 meV with the 
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FIG. SI: Temperature dependence of (circles) the specific heat 
Cp{T)/T of Bao.68Ko.32Fe2As2 and (green curve) its normal-state 
counterpart resulting from the corresponding state approximation by 
using the lattice specific heat of Ba(Feo.88Mno.i2)2As2 as a refer- 
ence, after Eq. (1) of the Letter 



FIG. S2: Temperature dependence of (circles) the specific heat 
Cp{T)/T of Bao.68Ko.32Fe2As2 and (blue curve) its normal-state 
counterpart resulting from the simulation of the lattice specific heat 
on a basis of Einstein modes, after Eq. (SI). Inset: The (5-functions 
Fi5{Lu — L0i) of the PDOS used for the fit are presented by a histogram. 



TABLE SI: Frequencies and weights of the six Einstein modes fit to the lattice specific heat of Bao.68Ko.32Fe2As2 (7 = 49 mJ/molK ). 



fkui , meV (K) 



1.92 
3.36 
5.89 
10.30 
18.04 
31.57 



(22.3) 

(39.0) 

(68.4) 

(119.6) 

(209.4) 

(366.4) 



0.005 
0.018 
0.136 
1.380 
1.450 
2.020 
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FIG. S3: Superconductvity-induced specific heat difference of of Bao.68Ko.32Fe2As2 obtained by subtractiongg the fit (green curve) by using 
the lattice specific heat of Ba(Feo.88Mno.i2)2As2 as a reference, after Eq. (1) of the Letter and Fig. Si, and (blue curve) by simulating the 
lattice specific heat on a basis of Einstein modes, after Eq. (SI) and Fig. S2. 

main contributions at about 10, 18 and 32 meV, that is in a very good agreement with measured by inelastic x-ray scattering 
and calculated phonon frequencies in doped and undoped double-layered BaFe2As2Compounds [S4]. The blue curve in Fig. S2 
represents the resulting normal-state counterpart of the specific heat of Bao.68Ko.32Fe2As2. Figure S3 compares results for the 
superconductivity-induced specific heat, ACei{T)/T, obtained as the difference between the measured specific heat Cp{T)/T 
of Bao.68Ko.32Fe2As2 and the normal-state specific heat resulting from two alternative approaches: by using the lattice specific 
heat of Ba(Feo.88Mno.i2)2As2 as a reference (open circles) or by simulating the lattice specific heat on a basis of Einstein modes 
(blue curve). The difference between these two curves are within the uncertainty represented by the symbol size. 

• Model spin-fluctuation spectral function, band and coupling parameters 

In our Letter, we desciibe the results of a fully microscopic calculation of the specific heat in the framework of a multiband 
spin-fluctuation model. The main input into the Eliashberg equations is the spectral function of the intermediate bosons. Early 
calculations of magnetic pairing in the high-Tc cuprate superconductors employed two alternative parametrized models of the 
susceptibility: The analysis of Radtke et al. [S5] [Radtke-Ullah-Levin-Norman (RULN) model] invoked neutron-scattering 
measurements, while that of Millis et al. [S6] [Millis-Monien-Pines (MMP) model] was based on NMR data. Because of a 
long tail at high frequencies cx l/w, the MMP spectrum leads to the divergence of the first momentum (proportional to the 
Hopfield-factor rj) and the overestimation of and superconducting gaps. Very recently, applying the MMP spectrum for spin 
fluctuations, the authors of Ref. S8 came to the conclusion about weak (intermediate) coupling in pnictides. The necessity of an 
improvement of the high energy dependence of the MMP spin-fluctuation spectral function was discussed in Ref. S7. In our 
analysis we took a more realistic coupling function based on neutron-scattering measurements, with a linear 51-dependence at low 
frequencies, a maximum at MVgf^ =18 meV, and a fast decay at > Ml™''^, in qualitative agreement with recent experimental 
data on the normal-state dynamical spin susceptibility of Ba(Fei_xCox)2As2 [S9, S 10]. With this spectral function shown in the 
insert of Fig.4 (also in Fig.S4) and the variable parameters described and listed in the Letter (also in Table SII), we obtained the 



TABLE SII: Microscopic paremeters of the four-band spin-fluctuation model with ftfjgj^'' = 18 meV and 13 meV fit the temperature dependence 
of the free-energy difference and yielding the corresponding Tc and gap A^'° values. N^'"^ is the density of states in the i*'' band, A is the 
four-band matrix of the coupling constants, and MSE is the mean square error of the fit. 
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FIG. S4: The best fit to the free-energy difference of the four-band Ehashberg model with the coupling function B{Q) peaked at hQ'l^p^ = 
18 meV (left) and 13 meV (right) and parameters listed in Table SII.. The red, green, and blue lines represent the partial contributions of the 
individual bands. 




FIG. S5: Temperature dependence of the superconducting energy gaps calculated from the four-band Eliashberg model with the coupling 
function B{^1) peaked at hfl'^p^ = 18 meV (left) and 13 meV (right) and parameters listed in Table SII. 



best simultaneous fit to the experimental values of Tc, the superconducting gaps, and the temperature dependence of the free- 
energy difference. The effective coupling constant averaged over all bands. A"" = J^rj K^''{0)Kj/Kot ~ 1-9' is remarkably 
close to the coupling constant of Pbo.8Bio.2 and confirms our conclusion about strong coupling in Bao.68Ko.32Fe2 As2 (see Fig. 3 
of the Letter for comparison). The consistency with the experiment tolerates some variation of the parameters and remains 
satisfactory for hil.'gp^ within the 10-20 meV energy range. In Table SII and Figs. S4-S5 we compare our results for two 
representative coupling functions with hilfp^ =18 meV and 13 meV, illustrating that a shift of hil'gp^ to lower energy leads to 
a larger values of A"" with corresponding densities of states iVj''''^(0) approaching the first-principle Density Functional results 
that accounts for the strong renormalization of the Sommerfeld constant 'y^/'joFT ^ 5. 
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